A formula for dielectric mixtures 
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Dielectric properties of material mixtures are of importance in diagnostics, characterization and 
design of systems in various engineering fields. In this Letter, we propose a peculiar dielectric 
mixture expression, which is based on the dielectric relaxation phenomena and the spectral density 
representation [E. Tuncer J. Phys. Condens. Matter 17(12) L125 (2005)]. The expression is tested 
on several composite systems. Results illustrate that the proposed expression can be used to obtain 
valuable structural informations in composites, even for highly filled, bi-percolating, systems. Lastly, 
the proposed expression is an alternative to other existing homogenization formulas in the literature. 

PACS numbers: 77.84.Lf 
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Predicting the dielectric properties of mixtures with 
or without knowing the composition of the mixture is an 
old problem in science and technology Q . The structure- 
property relationship in composite materials has been 
sought by researchers in various fields e.g., physics, elec- 
trical engineering, materials science, geophysics, pharma- 
ceutics, biology, etc., in order to comprehend the nature 
of material mixtures for many decades. Such a deep un- 
derstanding would be of great value to be able to cal- 
culate either the dielectric constant of a mixture of sub- 
stances of known dielectric constants or, knowing the di- 
electric constants of a mixture of two components and 
that of one of the components, to calculate the dielectric 
constant of the other [2| , or even knowing the dielectric 
constants of a mixture of two components and that of 
two components to estimate the morphology of the mix- 
ture Q • Recently, it has been shown Q that the applica- 
tion of the knowledge in dielectric relaxation phenomena 
and the theory developed by Bergman p| and Milton [(| , 
significant information could be obtained from the fre- 
quency dependent dielectric properties of a mixture when 
the dielectric properties of he constituents were known. 
In this Letter, we present a peculiar dielectric mixture 
formula, and utilize it to analyze the dielectric data 
of oil- in- water emulsion (O/W) Q, brine-porous rock 
(B/PR) @ and glass-beads-paraffin wax (GB/PW) @ 
systems. 

We first introduce a 'scaled' dielectric permittivity £ of 
the two-component mixture system as, 



£ = (£ e -£ m )(£ i -e m )" 



(1) 



where the subscripts 'e', 'm' and T denote the complex 
dielectric permittivities e with ohmic loss contributions 
of the effective, matrix and inclusion medium, respec- 
tively. Observe that £ varies between and 1. Second, 
the complex permittivities are expressed as, 

e k = e k x +x k M + a k (uve ) where k = {e, m,i} (2) 

Here and a are the dielectric permittivity at opti- 
cal frequencies and the ohmic conductivity of the ma- 
terial. The parameter x is the dielectric polarization 



30 
29 
3 s 28 



+ 



E 
CO 



27 

26 

25 
0.4 

0.3 

0.2 

0.1 









(a) 


oo 






o 


0000 oo-o— 












E 

a 1 

' a. 

w 05 














V 


(c) 






7* 











0.05 0.1 0.15 

9t(p [£J m]) 









9/ 

9^ 








* (b) 



















0.5 



1 1.5 
v [GHz] 



FIG. 1: Frequency dependent (a) relative dielectric permit- 
tivity and (b) conductivity a of the oil-in- water emulsion sys- 
tem. The Argand plot of the complex resistivity p is shown 
in (c). The measured and modeled data are presented with 
the open-symbols (o) and the solid line ( ), respectively. 



term, called the susceptibility, which is a complex and 
frequency v dependent quantity, ui in Eq. (J2J) is the angu- 
lar frequency (u = 2itv) , and i = \J— 1. We propose that 
the dielectric permittivity of a mixture can in principle 
be written as the dielectric relaxation with the 'scaled' 
permittivity replacing the dielectric permittivity of the 
material, Q 



i = 6 + A£[l + [wx) a ] 



ai-/3 



(3) 



'Electronic address: enis.tuncer@physics.org 



Here £ s is the fraction of inclusion phase 'i' forming a 
percolating network in the direction of the applied field — 
indicating the fraction of infinite cluster (s), which does 
not contribute to the interfacial polarization. A£ is the 
fraction of inclusion phase that is related to the iso- 
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TABLE I: The non-linear complex least-squares fit results. 
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O/W 


0.650 


3 x 10 -6 


4 x 10 -8 


0.962 


0.768 


B/PR 


0.064 


0.123 


0.032 
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0.741 


GB/PW 


0.025 


0.269 


0.152 


1 


1 



lated, deserted, clusters, which contribute to the polar- 
ization of the system. The permittivities of the con- 
stituents construct the so-called probing spectral fre- 
quency w(= e'/e m — 1). The quantity x is the spec- 
tral parameter, contains information on the effective de- 
polarization factors inclusions. The parameters a and 
f3 depend on the geometrical description of the system. 
Consequently the parameters x, a and f3 yield the dis- 
tribution of depolarization factors for the given system 
in hand|10l]. The concentration q of the inclusions are 
simply defined with q = £ s 4- A£; as a result for non- 
percolating systems q = A£ and completely percolating 
systems q = £ s . The defined parameters in Eq. Q are 
positive and smaller than one, < {£ s , A£, x, a, /?} < 1. 
Rearranging Eq. ©, we propose a dielectric mixing law 
as follows 



(e 1 - e m ){i s + A[l + {mxT]- } 



(4) 



We have selected three different systems to test the 
proposed expression in Eq. (@J; (i) oil- in- water 0, (it) 
brine-in-porous-rock and (Hi) metal-coated glass- 
beads-in-paraffin wax [9| composite systems. The data 
for the first two systems are digitized from the appropri- 
ate references. In the analysis we present the complex 
resistivity p[= (ieqEuj)' 1 ] level 0] of the systems as well 
to instruct the significance of this representation in lossy 
dielectric composites. 

Smith and Scott studied the dielectric properties 
of emulsion systems and applied the Bruggeman unsym- 
metric mixture formula [12] to investigate the behavior 
of their system, in which the oil was assumed to be 
the inclusion phase with lossless permittivity, and the 
frequency dependent complex permittivity of the water 
(saline solution) was calculated as presented in Ref. 
The permittivity of the composite system and response 
estimated with Eq. (gj by employing a non-linear com- 
plex least squares algorithm are shown in Fig.^ The fit 
parameters are listed in Table [I] The large discrepancy 
between measured and modeled data in Fig. QJi is due 
to the digitization procedure, however, the model values 
are within the range of the measured permittivity. The 
complex conductivity a estimate on the other is in ac- 
cord with the measurement and better than Bruggcman's 
expression^^ as shown in Ref. 0. Since the system shows 
conductive behavior, the Argand plot of the complex re- 
sistivity is also presented in Fig. to indicate the merit 
of the fitting expression. The concentration g(« 0.66) 
estimated indicates that the saline solution content is in- 
compatible with the value given by Smith and Scott 0, 
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FIG. 2: Frequency dependent (a) relative dielectric permit- 
tivity and (b) conductivity a of the brine-porous-rock system. 
The Argand plot of the complex resistivity p is shown in (c). 
The measured and modeled data are presented with the open- 
symbols (o) and the solid line ( ), respectively. 



q = 0.5. We therefore argue that although the mixture 
was composed of 1:1 oil and water, not all saline solu- 
tion was contributing to the conductivity as suggested 
by Smith and Scott [g . This observation means that the 
cmulsifier contant was higher than desired in the pre- 
pared mixture. 

Kenyon presented the dielectric data on White- 
stone saturated with salty water, the system was a brine- 
porous-rock mixture. Here we adapted the complex per- 
mittivities of the water and the rock taken from Stroud 
et al. 0. The fitting results are presented in Fig. [21 and 
the fit parameters are listed in Tabled It was stated that 
the porosity of the stone was 80% [= (1 — q) X 100] in Ref. 
H a value in accord with our finding w 0.813(= 1 — q = 
1 — £ s — A£) . In addition, even though the morphology of 
the rock was complex, the proposed expression gives evi- 
dence that some part of the brine was located in infinite 
pore clusters (£ s > 0). 

Youngs et al. |f| have recently presented the dielectric 
properties of paraffin wax filled with metal-coated hollow 
glass-beads on microwave frequencies. In the present in- 
vestigation, we adopt the complex permittivity of paraf- 
fin wax from the measurements, and estimate the inclu- 
sion permittivity from a dilute mixture measurement by 
utilizing an effective medium approximation |l4j . Ex- 
cept the data at low frequencies (y < 5 GHz), Eq. (gj 
yield results in accordance with the measurements. It 
is striking that the predicted concentration q(= 0.294) 
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FIG. 3: Frequency dependent (a) relative dielectric permit- 
tivity and (b) conductivity a of the glass-bead-praffin wax 
system. The Argand plot of the complex resistivity p is shown 
in (c). The measured and modeled data are presented with 
the open-symbols (o) and the solid line ( ), respectively. 

is in good agreement with the nominal volume fraction 
given by Youngs et al. The fit results listed in Ta- 
ble demostrate that most of the particles are arranged 



in an infinite cluster (£ s > A£) instead of forming isolated 
groups. 

In this Letter, we proposed a peculiar dielectric mix- 
ture expression, which do not suffer from the applicability 
and inaccurate predictions of previous mixing laws in the 
literature e.g. when inclusions form connected net- 
work, percolation, structure, composites with arbitrary 
shaped inclusions and highly filled mixtures. The ap- 
plication of the expression to several data sets has illus- 
trated its significance and potential in predicting and in 
comprehending the dielectric properties and structure of 
composites. We demonstrated that valuable information 
regarding the nature of mixtures can be obtained by the 
proposed expression. It can even be extended to systems 
exhibiting several depolarization processes, 

£ = & + £>^1 + (™*») att ] H \ ( 5 ) 

n 

which is analogous to theoretical representations of di- 
electric relaxation. Last but not least, it is suggested 
that the complex resistivity level is used while analyz- 
ing dielectric properties of composites jnj. Finally, the 
proposed expression can be applied to homogenization 
problems in other fields, e.g. elastic modulus, thermal 
conductivity, magnetic permeability, etc., where the ef- 
fective property of the heterogeneous medium is of im- 
portance. 

I would like to express my thanks to Drs. Nicola 
Bowler (Iowa State University, USA) and Ian J. Youngs 
(DSTL, UK) for supplying the dielectric data on the 
GB/PW system and for their fruitful comments. 
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